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Abstract 



We study a model of quasiparticles on a two-dimensional square lattice 
coupled to Gaussian distributed dynamical molecular fields. We consider two 
types of such fields, a vector molecular field that couples to the quasiparticle 
spin-density and a scalar field coupled to the quasiparticle number density. 
The model describes quasiparticles coupled to spin or charge fluctuations, 
and is solved by a Monte Carlo sampling of the molecular field distributions. 
When the molecular field correlations are sufficiently weak, the corrections to 
the self-consistent Eliashberg theory do not bring about qualitative changes 
in the quasiparticle spectrum. But for a range of model parameters near the 
magnetic boundary, we find that Migdal's theorem does not apply and the 
quasiparticle spectrum is qualitatively different from its Eliashberg approx- 
imation. In the range of model parameters studied, we find the transverse 
spin-fluctuation modes play a key role. While a pseudogap opens when quasi- 
particles are coupled to antiferromagnetic fluctuations, such a pseudogap is 
not observed in the corresponding charge-fluctuation case for the range of pa- 
rameters studied, where vertex corrections are found to effectively reduce the 
strength of the interaction. This suggests that one has to be closer to the bor- 
der of long-range order to observe pseudogap effects in the charge-fluctuation 



case than for a spin-fluctuation induced interaction under otherwise similar 
conditions. 

An important feature of the magnetic pseudogap found in the present 

calculations is that it is strongly anisotropic. It vanishes along the diagonal 

of the Brillouin zone and is large near the zone boundary. In the case of 

ferromagnetic fluctuations, we also find a range of model parameters with 

qualitative changes in the quasiparticle spectral function not captured by the 

one-loop approximation, in that the quasiparticle peak splits into two. We find 

that one needs to be closer to the magnetic boundary to observe the pseudogap 

effects in the nearly ferromagnetic case relative to the nearly antiferromagnetic 

one, under otherwise similar conditions. We provide intuitive arguments to 

explain the physical origin of the breakdown of Migdal's theorem. 
PACS Nos. 71.27.+a 
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I. INTRODUCTION 



The polarizer-analyser analogy provides an intuitive description of the effective inter- 
action between quasiparticles in a quantum many-body system. The first quasiparticle 
polarizes the medium in which it travels and the second quasiparticle, the analyser, feels the 
disturbance induced by the first one. In a strongly correlated system one can expect this 
induced polarization of the medium to be very complex, and in practice simplifying assump- 
tions are made. A commonly used approximation is the description of the polarization effects 
by the appropriate linear response function of the material. Furthermore, one typically only 
considers the interaction channel for which the linear response of the system is the largest. 
On the border of long-range magnetic order for instance, the spin-spin correlation function 
is the most enhanced and it is plausible that the dominant interaction channel is of magnetic 
origin and depends on the relative spin orientations of the interacting quasiparticles. 

It has been shown that such a magnetic interaction, treated in the self-consistent Eliash- 
berg approximation, can produce anomalous normal state properties and superconducting 
instabilities to anisotropic pairing states. It correctly predicted the symmetry of the Cooper 
state in the copper oxide superconductors [1] and is consistent with spin-triplet p-wave pair- 
ing in superfluid 3 He [for a recent review see, e.g., Ref. [2]]. One also gets the correct order of 
magnitude of the superconducting and superfluid transition temperature T c when the model 
parameters are inferred from experiments in the normal state of the above systems. For the 
case of a nearly half-filled single band, the calculations showed that the Eliashberg supercon- 
ducting transition temperature T c is higher for the tetragonal (quasi-2D) than cubic (3D) 
lattice [3-6]. Particularly striking is the comparison between the cubic antiferromagnetic 
metal Celn^ and the closely related compound CeCoIn 5 [7]. Superconductivity is found 
to extend over a much wider range in both temperature and pressure in CeCoIn^ than in 
Celnz- These findings and the growing evidence that the pairing symmetry in CeCoIn^ is 
d-wave in character, were correctly anticipated by the magnetic interaction model. 

While there have been a number of examples of superconductivity on the border of anti- 
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ferromagnetism, the corresponding phenomenon on the edge of metallic ferromagnetism has 
only been found recently. This result is not surprising within the framework of the magnetic 
interaction model. For otherwise equivalent conditions, the superconducting transition tem- 
perature is typically much higher on the border of antiferromagnetism than on the border 
of ferromagnetism [5,6,8]. An intuitive understanding of this finding pointed to candidate 
systems in which superconductivity on the border of ferromagnetism would be more likely 
to be observed. In particular, one ought to look for systems with strong spin anisotropy, 
i.e, with stong spin-orbit coupling and/or in a weakly spin-polarized state. This suggested 
a more detailed investigation of U Gt2 at high pressure which satisfied the above conditions 
and could be prepared in high purity form. This material proved to be the first example of 
the co-existence of superconductivity and itinerant-electron ferromagnetism [9]. 

The self-consistent Eliashberg treatment of the magnetic interaction model can produce 
strongly damped quasiparticles, but the electronic spectral function one obtains always shows 
a quasiparticle peak as one approaches the Fermi level. This is in qualitative disagreement 
with photoemission experiments on underdoped cuprate superconductors which show a near 
absence of a quasiparticle peak near the (vr, 0) point in the Brillouin zone [10]. This de- 
pletion of quasiparticle states, or pseudogap, is also seen in thermodynamic measurements 
[11]. The phenomenon may not be specific to the underdoped cuprates. For instance, the 
possible existence of a pseudogap in the heavy fermion compound CeCoIn 5 has recently 
been reported [12]. 

Is the disagreement between theory and experiment the manifestation of a fundamental 
flaw in the approach or does it simply reflect the inadequacy of the approximations used in 
the solution of the model? Questions regarding the validity of the Eliashberg treatment of the 
magnetic interaction model have been raised [13]. The one-loop approximation effectively 
assumes that quasiparticles behave as test particles. There must therefore be corrections to 
the simple theory, referred to as vertex corrections, coming from the fact that real and test 
particles behave differently. One can expect these vertex corrections to produce quantitative 
changes to the self-consistent Eliashberg theory as one approaches the border of long-range 
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magnetic order. In Ref. [14] it was shown that for optimally doped cuprates, vertex correc- 
tions did not bring about significant changes to the single spin-fluctuation approximation, 
producing an enhancement of the spin-fluctuation interaction of the order of 20%. The 
sign of this correction [14] is opposite to that expected in the case of a phonon-mediated 
interaction [15] and in the spin-density- wave phase [13]. This enhancement of the fermion 
spin-fluctuation vertex in the paramagnetic state is due to the transverse spin-fluctuation 
modes. 

It has been argued, in particular by Schrieffer [13], that as the antiferromagnetic corre- 
lations got stronger, the system should display characteristics akin to the antiferromagnetic 
insulating state and that the behavior of the quasiparticles should become qualitatively dif- 
ferent from that of a simple metal. The results presented here show that this physical insight 
is correct in that Midgal's theorem can qualitatively break down when the antiferromagnetic 
correlations become strong enough and a different state emerges. 

The validity of Baym-Kadanoff many-body theories such as the fluctuation exchange 
approximation [16], which is in many ways similar to the Eliashberg theory of the magnetic 
interaction model, has been extensively studied by Tremblay and collaborators in the con- 
text of the Hubbard model [17-19]. They find that close to the magnetic boundary, Migdal's 
theorem qualitatively breaks down, in that a critical-fluctuation-induced pseudogap (or pre- 
cursor pseudogap) is observed in the Quantum Monte Carlo simulations but is not found in 
the fluctuation exchange approximation. Similar results about the role of vertex corrections 
were reported in the case of the attractive Hubbard model [20-24], where in this case the 
precursor pseudogap is caused by critical pairing fluctuations. 

In this paper, we examine corrections to the single spin-fluctuation exchange approxi- 
mation in two dimensions using a non-perturbative formulation of the magnetic interaction 
model which is amenable to computer simulation [25] . We find that as the antiferromagnetic 
correlations become strong enough for vertex corrections to produce a qualitative change rel- 
ative to the one-loop approximation, a pseudogap opens in the quasiparticle spectrum. In 
this respect, our results are similar to those obtained for the Hubbard model [17-19]. The 
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pseudogap we find is strongly anisotropic in that it vanishes along the diagonal of the Bril- 
louin zone and is large near the zone boundary. We demonstrate that, in the range of model 
parameters studied here, the transverse spin-fluctuation modes are key to the appearance 
of the pseudogap by considering the case of commensurate charge-fluctuations with a spec- 
trum identical to that of the paramagnons. While a pseudogap opens when quasiparticles 
are coupled to magnetic fluctuations, such a pseudogap is not observed in the corresponding 
charge-fluctuation case for the range of parameters studied, where vertex corrections are 
found to effectively reduce the strength of the interaction. This suggests that one has to be 
closer to the border of long-range order to observe pseudogap effects in the charge-fluctuation 
case than for a spin-fluctuation induced interaction under otherwise similar conditions. In 
the case of nearly ferromagnetic systems, as magnetic correlations get stronger, we also find 
qualitative changes in the quasiparticle spectral function not captured by the one-loop ap- 
proximation. The quasiparticle peak splits into two distinct peaks with a lowering of the 
tunneling density of states at the Fermi level. For the range of parameters studied, this 
suppression of the tunneling density of states at the Fermi level is weaker than in the case 
of nearly antiferromagnetic systems for otherwise similar magnetic and electronic spectrum 
parameters. 

The paper is organized as follows. In the next section we describe the model and the class 
of vertex corrections considered. Section III contains the results of the numerical simulations. 
In section IV, we give intuitive arguments for the physical origin of the pseudogap. Finally 
we give a summary and outlook. Most of the technical details are included in an appendix. 

II. MODEL 

There are essentially two different ways to describe quantum many-body systems. In 
the Newtonian or Hamiltonian approach, one considers particles interacting with each other 
via pairwise interaction potentials. This is the point of view commonly adopted in the 
perturbation-theoretic approach to the non-relativistic many-electron problem. In the rela- 
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tivistic version of the theory, however, one adopts a Maxwellian point of view in which the 
interactions between electrons are mediated by a field, the quantized electromagnetic field. 
The Maxwellian approach is also widely used to carry out numerical simulations of inter- 
acting systems based on the Feynman path integral. In such functional integrals, bosonic 
fields are represented by c-numbers. But fermion fields must be represented by anticom- 
muting, or Grassmann, numbers that are not easily handled by digital computers. When 
a dynamical molecular (or Hubbard- Stratonovich) field is introduced to mediate the inter- 
actions between the fermions, the problem is reduced to that of non-interacting particles 
in a fluctuating field. The integrals over the anticommuting variables can then be evalu- 
ated exactly, at least formally, thus eliminating the troublesome Grassmann numbers from 
the problem. In the following, we adopt the Maxwellian point of view and use dynamical 
molecular fields to mediate the interactions between quasiparticles. 

In general, the distribution of Hubbard-Stratonovich fields is very complex and usually 
leads to intractable problems due to the infamous "fermion sign problem". In a many-body 
system, the presence of other particles produces changes in the effective interaction between 
two particles through screening effects and induces arbitrarily complex self-interactions of 
the dynamical molecular fields. 

In this paper, we assume that the renormalization of the effective two-body interaction 
can be accounted for by a redefinition of the parameters entering the bare interaction. And 
we ignore all the self-interactions of the Hubbard-Stratonovich fields. The model was intro- 
duced in Ref. [25] and bears some resemblance to the " quenched approximation" of lattice 
gauge theories introduced by Marinari et al., who, incidentally, chose the name by analogy 
to condensed matter physics [26]. A recent application of the quenched approximation to 
the pseudogap problem for static Hubbard-Stratonovich fields and a clear exposition of the 
formalism is given by Posazhennikova and Coleman [27]. But in the present work, however, 
we use dynamical rather than static Hubbard-Stratonovich fields and a non-separable form 
of the molecular field correlation function. 

Very close to the boundary of magnetic or charge long-range order, the self-interactions 
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of the dynamical molecular field ignored in the present work are known to be important for 
low dimensional systems [28]. The approximations made here may not be appropriate in 
a number of other cases. The virtue of the present approach is that it gives insights into 
quantum many-body problems that are essentially non-perturbative. More importantly, 
the results presented in this paper demonstrates that the simplest model already yields 
interesting physics. Some of the simplifications made here can in principle be relaxed and 
the theory extended accordingly. 

To be more specific, we consider particles on a two-dimensional square lattice whose 
Hamiltonian in the absence of interactions is 

ho(r) = -J2 tij^L(T)^ja(r) - l^^2^la(r)^ia(r) (2.1) 

i,j,<x %ol 

where tij is the tight-binding hopping matrix, fi the chemical potential and ijj} a , ipi a respec- 
tively create and annihilate a fermion of spin orientation a at site i. We take = t if sites 
i and j are nearest-neighbors and Uj = t' if sites i and j are next-nearest-neighbors. 

To introduce interactions between the particles, we couple them to a dynamical molecular 
(or Hubbard-Stratonovich) field. It is instructive to consider two different types of molecular 
fields. In the first instance, we consider a vector Hubbard-Stratonovich field that couples 
locally to the fermion spin density. This is the case considered in Ref. [25] and describes 
quasiparticles coupled to magnetic fluctuations. To illustrate the role played by transverse 
spin-fluctuations we shall also consider the case of a scalar field that couples locally to 
the fermion number density. This case corresponds to a coupling to charge-fluctuations or, 
within the approximation we are using here, "Ising"-like magnetic fluctuations where only 
longitudinal modes are present. The Hamiltonians at imaginary time r for particles coupled 
to the fluctuating exchange or scalar dynamical field are then 

h(r) = h (r) -^M^rj-iW^fr) (2.2) 
h(r) = ta-s^^toiiT) (2.3) 

ia 

where Mj(r) = (Mf(r), Mf (r), Mf(r)) T and $j(r) are the real vector exchange and scalar 
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Hubbard-Stratonovich fields respectively, and g the coupling constant. The reason for the 
choice of an extra factor l/\/3 in Eq. (2.2) will become clear later. 

Since we ignore the self-interactions of the molecular fields, their distribution is Gaussian 
and given by [25] 

V[M] = 1 exp ( - £ M(q,^) ; M(- q ,^ B )\ 

Z= /DMexpf-E^'^;^'"^! (2.5) 
J P V i& 2a(q,ii/ B ) ) 1 J 

in the case of a vector exchange molecular field and 

Z = f D$exp f - £ ^^M-q-^A (2.7) 

y p v q^x 2a(q,^z/ n ) ; 1 ; 

in the case of a scalar Hubbard-Stratonovich field. In both cases v n = 2imT since the 
dynamical molecular fields are periodic functions in the interval [0,(3 = 1/T]. The Fourier 
transforms of the molecular fields are defined as 

M r (t) = £ M(q, iv n ) exp ( - i[q • R - u n r]) (2.8) 
$r(t) = £ $(q, «/ n ) exp ( - i[q • R - i/ n r]) (2.9) 

We are considering the case where there is no long-range magnetic or charge order. The 
average of the dynamical molecular fields must then vanish and their Gaussian distributions 
Eqs. (2.4,2.6) are completely determined by their variance a(q, iv n ), which we take to be 

f iv n ) if M(q, iv n ) or $(q, iv n ) complex 

a(q> n ) = (2.10) 
I jfX(<l, w n) lf M (q> w n) or $(q, w «) real 

where iV is the number of allowed wavevectors in the Brillouin zone. Then 

T 

Mi(q,M/ n )Mj(k,in n )) = — x(q,w„)5 q -k^ n -n n ^j (2.11) 

T 

$(q, w/ n )$(k, ifi n )) = — x(q, Wn)<5q ) -k^„,-n„ (2.12) 
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where < . . . > denotes an average over the probability distributions Eq. (2.4) and Eq. (2.6) 
for the vector and scalar cases respectively. In order to compare the scalar and vector molec- 
ular fields, we take the same form for their correlation function x(q, iv n ) and parametrize it 
as in Refs. [5,8]. In what follows, we set the lattice spacing a to unity. For real frequencies, 
we have 

X(q,")= 2 J^. U (2-13) 

where k and k are the correlation wavevectors or inverse correlation lengths in units of the 
lattice spacing, with and without strong correlations, respectively. Let 

f ± = 4 ± 2(cos(g x ) + cos(g„)) (2.14) 

We shall consider commensuate charge-fluctuations and antiferromagnetic spin- 
fluctuations, in which case the parameters q 2 and rj(q) in Eq. (2.13) are defined as 

q 2 = q 2 + (2.15) 
v(q)=T q- (2.16) 

where T is a characteristic temperature. 

We shall also consider the case of ferromagnetic spin-fluctuations, where the parameters 
q 2 and 77(g) in Eq. (2.13) are given by 

q 2 = q 2 _ (2.17) 
v(q)=T q- (2.18) 

x(q, iu n ) is related to the imaginary part of the response function imx(q, a;), Eq. (2.13), via 
the spectral representation 

M)= .r-±!mM (2 . 19) 

J-00 7T %V n — UO 

To get x(q, iv n ) to decay as 1/v 2 as v n — > 00, as it should, we introduce a cutoff u and take 
Jmx(q, a;) = for u > u . A natural choice for the cutoff is uo = r)(q) Kq. 
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In the approximation we are considering, the single particle Green's function is the 
average over the probability distributions V[M] (Eq. (2.4)) or V[$\ (Eq. (2.6)) of the fermion 
Green's function in a dynamical vector or scalar field. 

g(iaT;ja'r') = J DMP[M] G(iar;ja'T'\[M\) (2.20) 

giiaT-jaW) = J D<$> V[§\ G(iar; 3<t't'\[§\) (2.21) 

where 

G(tar;ja'r'\[M] or [$]) = -(T r {^(r)^(r')}> (2.22) 

is the single particle Green's function in a dynamical molecular field and is discussed in 
the appendix. In evaluating expressions Eqs. (2.20,2.21) one is summing over all Feynman 
diagrams corresponding to spin or charge-fluctuation exchanges [25,27]. The diagrammatic 
expansion of the Green's function in a dynamical field, Eq. (2.22) and its average, Eq. (2.20) 
are shown pictorially in Fig. 1. A non-perturbative study of quasiparticles coupled to static 
magnetic fluctuations using a diagram summation technique was reported in Ref. [29]. In 
this paper we do not resort to a diagrammatic expansion but rather evaluate the averages 
in Eqs. (2.20,2.21) over the dynamical molecular fields by Monte Carlo sampling. 

It is very instructive to compare the results of the Monte Carlo simulations with the 
self-consistent Eliashberg calculations for the same model. If one only considers single spin 
or charge-fluctuation exchange processes, the single particle Green's function is given by: 

T 

E(p, iuj n ) = g 2 — Yl X(P - k > iuJ n - itt n )G(K itt n ) (2.23) 

S(P, iu n ) = ~. t \ w . ; (2.24) 

iu n - (e p - /I) - S(p, iu n ) 

where S(p, iu> n ) is the quasiparticle self-energy, Q(p,iu n ) the one-particle Green's function. 
e p is the tight-binding dispersion relation obtained from Fourier transforming the hopping 
matrix in Eq. (2.1) and jx the chemical potential. The choice of the factor 1 / \f?> in Eq. (2.2) 
means one obtains the same one-loop equations in the the case of an exchange molecular 
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field M and that of the scalar field $, thereby simplifying the comparison between the two 
cases. 

III. RESULTS 

The quasiparticle dispersion relation for the two-dimensional square lattice is obtained 
from Eq. (2.1). We measure all energies and temperatures in units of the nearest-neighbor 
hopping parameter t. We set the next-nearest-neighbor hopping parameter t' = — 0A5t. 
The chemical potential is adjusted so that the electronic band filling is n = 0.9. The 
dimensionless parameters describing the molecular field correlations are g 2 Xo/t, T /t, n Q 
and k. A complete exploration of the parameter space of the model is beyond the scope of 
this preliminary study. We chose a representative value for = 12, and set T = 0.67t as in 
our earlier work [5,8]. For an electronic bandwidth of leV, T ps 1000°K. We only consider 
one value of the coupling constant g 2 Xo/t — 2. In the random phase approximation, the 
magnetic instability would be obtained for a value of g 2 Xo/t of the order of 10. We consider 
what happens to the quasiparticle spectrum at a fixed temperature T = 0.25t as the inverse 
correlation length k changes. 

The calculations were done on a 16 by 16 lattice, with 41 imaginary time slices, or 
equivalently 41 Matsubara frequencies for the molecular fields, M(q, iv n ) and $(q, iu n ) {y n = 
2vmT, with n = 0, ±1, . . . , ±20). 

By analytic continuation of the single particle Green's function Q(k,r) or Q(k, iu n ) one 
can obtain the quasiparticle spectral function A(k,u) = —\lm Q R (k,u) and the tunneling 
density of states N(u) = X)k ^.(k, w), where Gr(\l, oj) is the retarded single particle Green's 
function. In the case of the one-loop approximation, Q(k, iuj n ) is analytically continued from 
imaginary to real frequencies by means of Pade approximants [30]. The imaginary time 
Monte Carlo data is analytically continued with the Maximum Entropy method [31]. We 
have used 2000 MC samples binned in groups of 20 to make 100 measurements. We have 
used three different versions of the Maximum Entropy method. The Classic MaxEnt and 
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two versions of Average MaxEnt, which is the method recommended in Ref. [31], where the 
probability for the a parameter is either a constant or proportional to 1/a [31]. In all cases 
we chose a flat default model. The results for the three different versions of the Maximum 
Entropy method tried were nearly identical in the case of the vector molecular field. There 
were slight differences between the Classic and Average MaxEnt solutions for the spectral 
function in the scalar molecular field case, while the two versions of the Average MaxEnt 
gave nearly identical results. The choice of the ~ 1/a probability distribution in the Average 
MaxEnt method gave a slightly better fit to the MC data and all the results shown here are 
those obtained with this choice of probability distribution for the a parameter. 

Fig. 2 shows the tunneling density of states N(u) and the spectral function A(k, oo) = 
— ^ImQ R (k,u) for the one-loop Eliashberg approximation to the Green's function, 
Eqs. (2.23,2.24), in the case of commensurate charge-fluctuations and antiferromagnetic 
spin-fluctuations. Figs. 2b, c show a strong anisotropy of the spectral function. A(k, uS) is 
sharper when k is along the diagonal compared to the case where k is near a hot spot 
(which is a point on the Fermi surface accessible from another via a momentum transfer of 
Q = (it, it)). Also note the monotonic broadening of A(k, u) as the parameter k 2 is reduced. 

In Fig. 3 we show the tunneling density of states and quasiparticle spectral function one 
obtains from the Monte Carlo calculation with a coupling to the scalar dynamical molecular 
field $. By comparison to the one-loop self-consistent results of Fig. 2, the vertex corrections 
give rise to a sharpening of the quasiparticle spectral function, except for k 2 = 4 and 
k = (37r/8, 37r/8). For the range of values of model parameters considered here, the multiple 
charge-fluctuation exchanges do not lead to a breakdown of the quasiparticle picture. Note 
that contrary to the Eliashberg result, as k 2 is reduced from 4, the spectral function initially 
sharpens before broadening again at the lower values of k 2 considered here. For k 2 = 4, 2, 
the spectral function is sharper near the hot spot than along the diagonal, again in contrast 
to the one-loop self-consistent result. It is approximately isotropic at k 2 = 1 and becomes 
sharper along the diagonal than near the hot spot at the lower values of k 2 = 0.50, 0.25. 
Quite generally Fig. 3 shows that when vertex corrections are included the spectral function 
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is less anisotropic than the Eliashberg result for the values of k 2 considered. 

The Monte Carlo results for quasiparticles coupled to antiferromagnetic spin-fluctuations 
are shown in Fig. 4. For k 2 < 1, a pseudogap appears at the point k = (jv,n/8) but not 
along the diagonal, in qualitative agreement with experiments on underdoped cuprates [10]. 
The pseudogap also shows up in the tunneling density of states N(cu) which is suppressed at 
the Fermi level for k 2 < 1. By comparing the spectral function of Figs. 4b, c to the one-loop 
self-consistent result, we see that for the values of k 2 where there isn't a pseudogap, vertex 
corrections lead to a broadening of the spectral function and a reduction of the momentum 
anisotropy of the spectral function. The contrast between the non-perturbative calculations 
for charge and spin-fluctuations and the comparison with the Eliashberg result is shown in 
Fig. 5 for k 2 = 2. In this case it is clear that the effect of vertex corrections is to make 
the spectral function sharper than its Eliashberg approximation in the charge-fluctuation 
case and broader when the quasiparticles are coupled to magnetic fluctuations. For these 
model parameters, corrections to the Eliashberg approximation reduce the effective charge- 
mediated interaction and enhance the coupling to magnetic fluctuations. 

In Fig. 6 the results of the one-loop calculations in the case of ferromagnetic spin- 
fluctuation exchange are shown. As could be expected, the spectral function for quasi- 
particles near the Fermi level does not depend strongly on whether the momentum is near 
the hot spot, Fig. 6b, or along the diagonal in the Brillouin zone, Fig. 6c, in contrast to the 
nearly antiferromagnetic case. At the Eliashberg level, the spectral function broadens mono- 
tonically as the correlation wavevector k is reduced and is broader than the corresponding 
nearly antiferromagnetic case, Fig. 2b, c. 

The Monte Carlo results for quasiparticles coupled to ferromagnetic spin-fluctuations 
are shown in Fig. 7. The tunneling density of states at the Fermi level N(uj = 0) begins 
to drop as the parameter k 2 < 0.50. For k 2 = 0.25, the quasiparticle peak has effectively 
split into two peaks, and a precursor of this effect can be seen at k 2 = 0.50. This phe- 
nomenon is qualitatively different from what one can obtain at the one-loop level, Fig. 6, 
just as the pseudogap seen in Fig. 4, when the quasiparticle interactions are mediated by 
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an antiferromagnetically correlated exchange field. For the larger value of k 2 where the 
splitting of the quasiparticle peak is not observed, a comparison of the Monte Carlo results 
with the one-loop self-consistent calculations shown in Fig. 6 shows that vertex corrections 
bring about a broadening of the quasiparticle spectral function. Therefore, just as in the 
antiferromagnetic case, vertex corrections enhance the magnetically-mediated interaction, 
albeit to a lesser degree, as a close look at Figs. 2,4 and 6,7 indicates. This could explain 
why qualitative changes from the Eliashberg solution are seen for smaller values of k 2 for 
coupling to ferromagnetic than antiferromagnetic fluctuations. 

IV. DISCUSSION 

For the model considered here, at the one-loop level the exchange of commensurate 
charge-fluctuations and antiferromagnetic spin-fluctuations yield the same quasiparticle 
properties. Once vertex corrections are included, the differences between the two cases 
are evident. In all but one case, k = (37r/8, 37r/8) and k 2 = 4, corrections to the one-loop 
approximation make the quasiparticle peak sharper than the Eliashberg result for charge- 
fluctuations. In the case of antiferromagnetic spin-fluctuations, the quasiparticle peak gets 
broader as vertex corrections are included and a pseudogap appears for the lower values of k 2 
studied. The difference between the scalar and vector molecular fields appears at the leading 
vertex correction to the one-loop approximation [14] shown in Fig. 8. The frequency and 
momentum integrals are the same in both cases since we assume the commensurate charge- 
fluctuations and antiferromagnetic spin-fluctuations have the same spectrum. However, the 
spin sums in the two cases are not identical. In the case of a coupling of the molecular 
field to the quasiparticle spin-density, one gets a factor coming from the Pauli matrices at 
each vertex in the diagram. For the diagram shown in Fig. 8, this factor is J2i,j o-' l (T j cr l a j . 
One can split the sum into the i = j and i ^ j terms, and use the fact that a l a l = 1 and 
crV = -crV if % ^ j. Then £ij crVoV = E^VVV - £^ trVVV = 3 - 6 = -3. 
Note that the longitudinal spin-fluctuations contribute a term a z a z a z a z = 1 and thus the 
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change in sign is caused by the presence of transverse magnetic modes. The corresponding 
factor in the case of charge- fluctuations is 1, just as in the case without transverse mag- 
netic modes. Coupling the quasiparticles to the spin-density instead of the number density 
produces a leading vertex correction with the opposite sign. The process depicted in Fig. 8 
enhances the magnetic interaction [14] while this same diagram leads to a suppression of the 
effective interaction in the case of charge-fluctuations. For the range of parameters stud- 
ied, the Monte Carlo simulations show that this qualitative difference between coupling to 
magnetic or charge fluctuations persists to all orders, save one case, k = (37r/8, 37r/8) and 
k 2 = 4, where vertex corrections enhance the charge-mediated interaction. 

As noted above, for large values of k 2 , vertex corrections do not produce qualitative 
changes, but merely quantitative ones. One might therefore ask whether one can obtain the 
non-perturbative results with a one-loop calculation provided the parameters of the theory 
are renormalized. To illustrate the point, let us focus on the case of antiferromagnetic 
paramagnons with k 2 = 2. One essentially has two parameters one can renormalize, the 
dimensionless coupling constant g 2 Xo/t an d the inverse correlation length k. Since vertex 
corrections make the spectral function broader and the quasiparticle residue smaller, one 
could attempt to fit the non-perturbative calculations with a one-loop theory with a larger 
coupling constant g 2 Xo/t or smaller k 2 , or a combination of both. If the corrections to 
the one-loop theory are purely local, it is possible to absorb them in a redefinition of the 
coupling constant g 2 Xo/t- It turns out that increasing the coupling constant does make the 
quasiparticle residue smaller but the one-loop spectral function remains too sharp relative 
to the non-perturbative calculation. By the time g 2 Xo/t is large enough for the quasiparticle 
lifetime to be approximately that obtained with the Monte Carlo simulation, the one-loop 
quasiparticle residue is then too small. One has more success with making k 2 smaller and 
Fig. 9 compares the spectral function of the one-loop calculation with k 2 = 0.25 with the 
A(k, uo) obtained from the Monte Carlo simulations at k 2 = 2. One can get a decent fit near 
the hot spot at k = (tt, 7r/8) but the one-loop spectral function is always more anisotropic 
in momentum than the non-perturbative A(k,uj) for this value of k 2 . 
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The most important feature of the calculations presented here is the qualitative change 
in the quasiparticle spectral function that occurs as the magnetic interaction gets stronger, 
either on the border of long-range antiferromagnetic or ferromagnetic order. The appearance 
of a pseudogap in the quasiparticle spectrum near a second-order phase transition as the 
temperature approaches the critical transition temperature from above was demonstrated 
for the half-filled Hubbard model [17-19], for a superconducting instability [18,20-23,33,34] 
and for a Peierls-CDW transition [32]. By contrast, in the model studied in this paper, the 
quadratic actions of the dynamical molecular fields do not exhibit a phase transition at all 
(with k 2 > 0). It may therefore be somewhat surprising that a pseudogap is observed in 
the present calculations as the correlation wavevector k becomes of order one in the nearly 
antiferromagnetic case and of order one-half for nearly ferromagnetic systems. 

The physical origin of the pseudogap was explained in Refs. [17-19] and later in Ref. [32]. 
Quasiparticles only remain coherent for a finite amount of time. When the distance they 
can travel during that time becomes shorter than the correlation length of the molecular 
field, quasiparticles effectively see long-range order. In the presence of long-range anti- 
ferromagnetic order, the spin-density-wave quasiparticle spectral function consists of two 
peaks. In the case of long-range ferromagnetic order the spin-up quasiparticles have an 
energy shifted downwards, say, relative to the paramagnetic quasiparticle energy while the 
energy of a spin-down quasiparticle is shifted upwards relative to its energy in the absence of 
long-range ferromagnetic order. The shifts in energy give rise to corresponding shifts in the 
quasiparticle peak in the spectral function. Upon averaging over the two spin orientations, 
the quasiparticle spectral function would then also consist of two peaks, the up-spin and 
down-spin quasiparticle peaks. In our calculations, when the quasiparticles remain coherent 
for such a short time that they effectively see long-range order, the spectral functions have 
charactersitics akin to that of the ordered state. For instance, in Fig. 4b, for a coupling 
of quasiparticles to antiferromagnetic spin-fluctuations, the spectral function A(k, uj) for 
k 2 = 0.25 looks like that of a broadened spin-density wave. And in the case of a ferromag- 
netically correlated molecular field, in Figs. 7b, c, A(k,u) is "spin-split" as k 2 is reduced, a 
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feature which can be understood if, during their short lifetime, the quasiparticles effectively 
see ferromagnetic order, where the moment is equally likely to point up or down, since we 
are still in the paramagnetic phase. 

Quasiparticles only remain coherent for a finite time due to thermal and quantum fluctu- 
ations. The question is what is that time scale, or the associated characteristic length scale 
to be compared against the magnetic correlation length. In the half-filled Hubbard model 
studied in Refs. [17-19], the renormalized classical regime for the spin-fluctuations always 
precedes the zero temperature phase transition, and in that regime thermal fluctuations 
dominate and the relevant charactersitic length scale is the thermal quasiparticle de Broglie 
wavelength £ th = v F /T [17-19]. Thermal fluctuations are dominant near the Peierls-CDW 
transition in two dimensions and the relevant quasiparticle length scale for the onset of the 
pseudogap also turns out to be £ th [32]. For the model considered in this paper, however, 
thermal fluctuations do not appear to dominate. Indeed, the onset of the pseudogap for 
coupling of quasiparticles to antiferromagnetic fluctuations and the qualitative change seen 
in the spectral function for coupling to ferromagnetic fluctuations occur for £ th K > 1. The 
departure from the criterion £ th K ~ 1 for the breakdown of the Midgal approximation is 
larger in the ferromagnetic case. The above suggests that another, shorter, length scale is 
relevant in the present case. A possible candidate length scale is the quasiparticle mean-free 
path. 

One can attempt to make this more quantitative by extracting quasiparticle lifetimes 
and mean-free paths from our numerical results. In order to do this we fit the quasiparticle 
peak of the spectral function with a Lorentzian A L (k,u) = ^ (ui-E^+r' 2 w ^ich corresponds 
to a quasiparticle approximation for the retarded Green's function Qnik, cu) = CJ _ £ f k _ iFk 
describing the propagation of quasiparticles of energy E^. T k is related to the quasiparticle 
lifetime r k through T k = If one ignores the momentum dependence of the self-energy, 
the quasiparticle residue z k is related to the ratio of the band to the effective mass z k 
^pr. We define the bare velocity = ^(ff 1 ) 2 + (ft 1 ) 2 where e k is the band dispersion 
relation and the renormalized velocity v k = Zkvj® . The quasiparticle mean- free path is then 
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approximately Zk = v^r^. The values of Zk one obtains for charge-fluctuations as well as 
antiferromagnetic and ferromagnetic spin-fluctuations are shown in Fig. 10. In the case of 
a coupling of quasiparticles to charge-fluctuations, Fig. 10a, Zk for all the values of k 

considered. In this case, the quasiparticles travel far enough during their lifetime to see there 
is no long-range charge order. In the case of antiferromagnetic spin-fluctuations, Fig. 10b, 
one sees that at k = (tt,7t/8) Zk ~ 1/k at k 2 = 1, and that is where the pseudogap begins 
to appear in the spectral function, Fig. 4b. At momentum k = (37r/8, 37r/8), the mean-free 
path Zk > 1/^ at k 2 = 1 and one wouldn't expect a pseudogap. Zk becomes less than the 
correlation length at the lowest value of k 2 and one can see hints of a developing pseudogap 
in the spectral function, Fig. 4c. In the ferromagnetic case, Fig. 10c, the appearance of the 
two peaks in the spectral function is broadly consistent with the condition Z k < 1/k. Note 
that we are using a simple criterion, in an attempt to capture the essential aspects of the 
problem, to understand the emergence of a pseudogap in our calculations. The crossover to 
the new state is likely to depend on other details not taken into account by our criterion and 
is therefore not expected to occur exactly at Z k = 1/k. Moreover, the tails in most of the 
spectral functions are not Lorentzians and thus our definition of the quasiparticle lifetime is 
clearly approximate. 

We have also extracted mean-free paths for the one-loop approximation using the same 
methodology (results not shown). For low enough values of k 2 , one can get in the regime 
where Zk < 1/k, but no pseudogap in the spectral function is observed. What is the one- 
loop approximation missing? A candidate explanation, based on an analogy, is the follow- 
ing. When treating the potential scattering of a particle in a momentum space basis, to 
lowest order of perturbation theory (Born approximation), one assumes the wavefunction is 
unchanged, i.e remains a plane-wave. In order to study bound states, in which the wavefunc- 
tion of the particle is qualitatively different since it is localized, one must treat the scattering 
events to all orders. By analogy, the Eliashberg approximation does not seem to allow for 
a change in the quasiparticle wavefunctions which essentially remain plane-waves. In the 
pseudogap state, the quasiparticle wavefunctions must be qualitiatively different and one 
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must allow the quasiparticles to scatter multiple times against the locally ordered molecular 
field to produce the required changes in their wavef unctions. It seems one would need to 
sum an infinite set of spin-fluctuation exchanges. 

The "quenched approximation" formulation of the magnetic interaction model actually 
goes beyond a diagrammatic perturbation expansion. The Green's function, Eq. (2.20), 
depends on the position, spin, and imaginary time, but also on the parameters of the theory, 
A = g 2 Xo/t, k 2 , etc., Q = Q(iaT;ja'r'\\, k 2 , . . .). (To make the dependence on A = g 2 Xo/t 
explicit, one simply needs to introduce scaled molecular fields m = gM.. The variance of 
the new variables is then g 2 a(q,iu n ) oc # 2 x(q, iv n )). It is clear from Eqs. (2.20,2.4,2.5) 
that the theory does not make sense for A < 0, since the Gaussian distributions of the 
molecular fields would have a negative variance. This observation means that it is not 
possible to analytically continue Q to negative values of A, implying an essential singularity 
at A = 0. The perturbation expansion in powers of A = is then an asymptotic rather 
than convergent series [35]. This opens the possibility for phenomena that lie outside of 
diagrammatic perturbation theory. Whether the pseudogap state found in the numerical 
simulations reported on in this paper is precisely one such phenomenon is not presently 
known to the author. 

The results presented here do not imply one couldn't get a pseudogap state when coupling 
to charge fluctuations. On the basis of the arguments presented above, if one were to increase 
the strength of the charge correlations or the coupling of quasiparticles to the molecular field 
$, one would get in the regime C 1/k and the spectral function ought to resemble that 
of a broadened charge-density-wave state. Our calculations simply show that one would 
have to be closer to the ordered state in the case of charge fluctuations than for magnetic 
fluctuations, under otherwise similar conditions. 
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V. OUTLOOK 



We studied a non-perturbative formulation of the magnetic interaction model, in which 
quasiparticles are coupled to a Gaussian distributed dynamical molecular exchange field. 
Far from the magnetic boundary, the type of vertex corrections considered here do not 
bring about qualitative changes to the quasiparticle spectrum. But as one gets closer to 
the border of long-range magnetic order, we find, for a range of model parameters, that 
Migdal's theorem does not apply and the quasiparticle spectrum is qualitatively different 
from its Eliashberg approximation. The physical origin of the phenomenon is that if the 
distance quasiparticles can travel during their lifetime becomes shorter than the molecular 
field correlation length, these quasiparticles effectively see long-range order. When the 
molecular field correlations are antiferromagnetic, the quasiparticle spectral function has 
the two-peak structure of a spin-density-wave state, even though there is no spontaneous 
symmetry breaking. We find that the associated pseudogap is strongly anisotropic in that 
it vanishes along the diagonal of the Brillouin zone and is large near the zone boundary, in 
qualitative agreement with photoemission experiments on underdoped cuprates [10]. The 
anisotropy of the pseudogap found in our calculations simply reflects the anisotropy of the 
quasiparticle mean-free path. For coupling to ferromagnetic fluctuations, we also find a 
range of parameters were the quasiparticle spectral function becomes qualitatively different 
from its one-loop self-consistent approximation. The local ferromagnetic order also leads to 
a splitting of the quasiparticle peak into two. These pseudogap effects are found to be weaker 
for nearly ferromagnetic systems than for their nearly antiferromagnetic counterparts, under 
otherwise similar conditions. 

In the standard theory of quantum critical phenomena [28], if d + z < 4 where d is the 
spatial dimensionality and z the dynamical exponent, the mode-mode coupling parameter 
diverges upon renormalization as one approaches the instability and the critical exponents 
are different from their mean-field values. Therefore, if d + z < 4, one would certainly 
expect Migdal's theorem to qualitatively break down close enough to the quantum critical 
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point. However, when d + z > 4, the critical exponents take their mean-field values, as in a 
one- loop calculation. In that case, one could therefore expect, and it is often assumed, that 
the Eliashberg approximation is at least qualitatively correct. For the magnetic spectrum 
considered here, Eq. (2.13), z = 2 for antiferromagnetic fluctuations and z = 3 for ferromag- 
netic fluctuations. Hence d + z — 4, the marginal dimension in the case of antiferromagnetic 
spin-fluctuations but crucially, d + z = 5 in the case of ferromagnetic spin-fluctuations. One 
could thus have expected that in the latter case, the one-loop self-consistent calculations 
should be at least qualitatively correct since d + z > 4. This is at variance with our re- 
sults, which show qualitative differences bewteen the non-perturbative calculations and the 
Eliashberg predictions. Note that the qualitative break down of Migdal's theorem occurs for 
smaller value of n 2 in the ferromagnetic case than for antiferromagnetic spin-fluctuations. 
In that sense the effect of vertex corrections is weaker when d + z = 5 compared to the case 
d + z — 4, which is what is expected. 

To summarize, an often assumed criterion for the qualitative applicability of the Eliash- 
berg theory near an instability, namely d + z > 4, is clearly a necessary condition but the 
calculations presented here show that it is not a sufficient one. This result may point to 
certain limitations of the standard theory of quantum critical phenomena [28], which has 
recently been criticized by Anderson [36]. 

The crucial role of dimensionality for pseudogap phenomena has been emphasized by 
Tremblay and collaborators [17,18,21] and by Preosti et al. [37]. Since critical fluctuations 
responsible for the precursor pseudogap are much stronger in two dimensions than for three 
dimensional systems, they find that pseudogap phenomena are much weaker in three dimen- 
sions. One would expect similar results for the magnetic interaction model studied in the 
present paper. The role of lattice anisotropy is already important at the Eliashberg level 
[3-6], because of the increased phase space of soft magnetic fluctuations in lower dimensional 
systems. This phase space argument leads one to expect weaker vertex corrections to the 
Eliashberg theory of the magnetic interaction model, and hence weaker pseudogap effects in 
isotropic three dimensional systems. 
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The results presented here raise a number of obvious questions, which we hope to answer 
in the future. For instance, what are the corrections to the Eliashberg theory for the super- 
conducting transition temperature to the d x 2_ y 2 pairing state? We have studied corrections 
to quasiparticle spectral properties, and found that vertex corrections can produce a pseudo- 
gap in the quasiparticle spectrum. When the magnetic correlations are weak enough vertex 
corrections do not bring about qualitative changes but nevertheless reduce the quasiparticle 
lifetime. The above effects are expected to suppress magnetic pairing. But one must also 
include in the calculation of the superconducting transition temperature the corresponding 
corrections to the pairing interaction. In leading order, it was shown [14] that vertex cor- 
rections lead to a stronger pairing interaction in the d-wave channel. There may therefore 
be some cancellation of errors, at least in some range of model parameters. At this stage, 
however, one can only speculate about the effect of vertex corrections on the Eliashberg 
theory of the superconducting transition temperature. 

Transport and thermodynamic properties in the quenched approximation would also be of 
great interest. The calculations in Ref. [32] showed that the pseudogap in the single particle 
density of states also appeared in the two-particle Green's functions such as the optical 
conductivity and the uniform Pauli susceptibility. I therefore expect that the pseudogap 
observed in the quasiparticle spectral function will show up in the thermodynamic and 
transport properties of the model studied here. The extent to which the model is able to 
explain the normal state experimental data on the underdoped cuprates is an open question. 

Our simplifying assumptions, which were born out of the necessity to carry out the cal- 
culations in a reasonable amount of time, should be relaxed. The dynamical molecular field 
correlation function which enters their distribution was taken to be of the same functional 
form for all the calculations. In Nature, however, the distribution function of the molecu- 
lar field experienced by one quasiparticle is determined self-consistently by all of the other 
quasiparticles in the system. It is therefore expected to change as the nature of the quasi- 
particle spectrum changes. One should also study the effect of mode-mode coupling terms 
which were ignored in the present study. Monien has recently carried out such a study for 
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a model of the one-dimensional Peierls-CDW [38]. 

There is no doubt that the extension of the theory to deal with effects ignored here will 
bring about quantitative changes to our results. It will be of great interest to figure out if 
and for what model parameters our results are modified qualitatively. I would like to think 
that the physical origin of the emergence of a pseudogap in the quasiparticle spectrum will 
turn out to be independent of the details of the model. 
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VII. APPENDIX: GREEN'S FUNCTIONS 

In this appendix we derive the mathematical formulas for the single-particle Green's 
functions in a fluctuating exchange or scalar field. We also give the algorithm for the 
numerically stable and efficient calculation of such Green's functions. We shall make use 
of fermion coherent states and will derive the path-integral for an anti-normal ordered [39] 
Hamiltonian quadratic in the fermion field operators. 

We first need an identity for expressing the exponential of the quadratic form 

in terms of fermion coherent states. In the above expression, A = A^ is a Hermitian matrix 
of dimension N and {if)}, {ip } are fermion operators. One has: 

r N 

exp ( £ iM^J) = J [II C^i] le > exp ( E < £1 (7-2) 

where {£*}, {£} are Grassmann numbers. This identity is proved by going into a basis in 
which the matrix A is diagonal. Let U be the unitary transformation that diagonalizes A 
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rfAU = A = diag(Ai, . . . , X N ) (7.3) 
In terms of the new fermion fields {$}, {&} 



h = E utj*s ( 7 - 4 ) 

3 

4 = ]T Uu*] (7.5) 

3 

** = E ^ (7-6) 

3 

*1 = E (7-7) 

3 

the quadratic form Eq. (7.1) becomes Q A = EiAi$i$j. Since ($«$•)" = for n = 

1, 2, . . . , oo, and the terms with different indices i commute with one another, one has 

exp(Q A ) = JJeM^M) = 1] {l + [exp(A.) - 1]$,$!} (7.8) 

i i 

Expanding the product into sums yields 

exp(Q A ) = 1 + £[exp(AiJ _ 1]$^ + £ [exp(A il ) - l][exp(A i2 ) - 1]^^^^ 

il ii<i-2 

+ J2 i e MK) - l][exp(A i2 ) - l][exp(A l3 ) - l^^^iA^a + . . . 

h<i2<i3 

+ £ [exp(A il )-l]...[exp(A i J-l]^ 1 <...$ ijv $J JV (7.9) 

ii<...<ijv 

Since in the various sums one has i\ ^ i 2 ^ H ■ ■ ■ 7^ ^n, one can straightforwardly anti-normal 
order the terms 

^M^M 2 = ^MM, (7.io) 

$. $t d>. cbt d>. <|>t = $. $. $ cbt cbt cbt f7 in 

$, $ t $, <3>t $ f = $„ $, $ t $ f $t f7 12) 

We can now insert a resolution of the identity between the annihilation and creation opera- 
tors. For instance 

r N 

q^^MMMs = ****** /[nc<«>exp <^i«< 

J i=l i 
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/ nc^a^Cuie > exp ( - < m*c 2 c 3 

J i=l i 
N 

tn c^i 6,^*^* exp ( - x: ie >< ei 

i=l i 

Applying the above operation to every term in the sum in Eq. (7.9) one obtains 

r N N 

e W (Q A ) = j > exp ( lit 1 + ( eAl " < £1 

i=i i i=i 

. v 

= / [II dfrfc] le > exp ( - £ exp ( £(e A * - < £1 

j=i i i 

. v 

= / [II > exp (£e A ^*) < CI (7-13) 

i=l i 

Going back to the original representation in which the matrix A is not diagonal yields 
Eq. (7.2). If we have another quadratic form 

q b = ^M) ( 7 - 14 ) 

where B = is a Hermitian matrix of dimension N, it is a simple exercise in Grassmann 
integration to show that 

r N 

exp(Q A ) exp(Q B ) = / [JJ > exp ( ]T &( e A e s )^*) < CI (7.15) 

i=l ij 

To prove the above identity one uses the coherent state representation Eq. (7.2) for exp(<5yi) 
and exp(Q B ) 



r N 

exp(Q A ) = / [Ud^dm > exp (£&( e %$) < CI 

i=l ij 



iV 



(7.16) 
(7.17) 



exp(Qs) = y [n^*^]^ > ex P (^2vi(e B )i,jVj) < v\ 

i=l 

Integrating over the Grassmann variables {£*} and {r?;} yields Eq. (7.15). 

The imaginary time single-particle Green's function in a time dependent field is given by 
the usual expression 



G(tar;ja'r') = -(T r {^ CT (r) <,(r')}> 



Z = Tr 



-±Tr{T T { exp ( - J* h(s)ds)^ a (r)^ ja/ (r f )} 
T T exp ( — y h(s)ds^) 



(7.18) 
(7.19) 
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where the time dependent Hamiltonian h(r) is given by Eq. (2.2) and Eq. (2.3) for coupling 
to a fluctuating exchange or scalar field respectively. The trace is over the fermion fields. For 
ease of notation, we have dropped the explicit dependence of G(iar; ja'r') on the molecular 
field. Because of the time ordering operator, one must distinguish the two cases r > r' and 
t < t'. 



G(iar;ja'r') = ~^Tr 



{T T exp ( - h(s)ds)}iP ia (r){T T exp ( - jT h(s)ds)}^] a ,(T r ) 



x {T r exp ( - £ h(s)ds)} 



T > T 



(7.20) 



G(mr; joV ) = -^Tr 



{T r exp(- £h(s)ds)}^l,(r'){T T exp(- jT h(s)ds)}Aa(r) 



x {T r exp(-y o h(s)ds)} 



t <r' 



(7.21) 



The imaginary times r and r' are discretized such that r = Arm m = 0,1,... L t , 
t' = Arm' m' = 0,1,... L t , where At = One then uses the Trotter break- 

up for the Hamiltonian in an exchange, Eq. (2.2), or scalar field, Eq. (2.3), and hence 
exp(-Ar/f(mAr)) w exp(— Ar^i(mAr)) exp(— Ar^o(mAr)). Up to an unimportant con- 
stant, the anti-normal ordered hopping Hamiltonian ho, Eq. (2.1) is written as 



Tier J a' {tij ~\~ ^$i,j)$o,o' 



(7.22) 
(7.23) 



Similarly, up to constant terms that drop out of the expression for the single particle Green's 
function, the Hamiltonian for the coupling to the fluctuating field is written as 



(7.24) 



la j a' 
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A i<T , ja ,(m) = ^^[M^mAT)^^,! - M*{mAT)5„ A 5 a , A 

+ (Mf(mAr) + iM?(mAT))5 aA 5 a > A + (Mf(mAr) - zM^mAr))^^] (7.25) 
A i<Tda ,(m) = gSijS^^iimAr) (7.26) 

where Eq. (7.25) and Eq. (7.26) are for a fluctuating exchange and scalar field respectively. 
Making use of Eqs. (7.2,7.15), one has 
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exp(— Ar^(mAr)) « exp(— Ar/ii(mAr)) exp(— Ar/i (mAr)) 

= / ^(OIC > exp ( £ U^mk^) < CI (7.27) 

£ m = exp(-ArA(m)) exp(-ArT) (7.28) 

where from now on we use d/j,(£) to denote the Grassmann integration measure 

^(^[JICA] (7.29) 
jo- 
in the case of coupling to a scalar field 

[exp(-ArA(m))] iff)J v = <^<W exp(-#Ar$i(mAr)) (7.30) 
The case of the exchange field is slightly more complicated: 



[exp(-ATA(m))] iaJa , = 


<J i j[exp(-Ara i (m))] <T>ff / 


(7.31) 


exp(— Araj(m)) = 


/c - sA4f (mAr) -s.M+(mAr) \ 
\ -s.M~(mAr) c + sMfimAr)) 


(7.32) 


c = 


cosh^ Ar | Mi (mAr) | / v 7 ^) 


(7.33) 


s = 


sinh^Ar | Mi (mAr) | /v 7 ^) 


(7.34) 


M-(mAr) = 


Mf(mAr) 
|Mi(mAr)| 


(7.35) 


Mf(mAr) = 


Mf(mAr) ±iMf(mAr) 
|Mi(mAr)| 


(7.36) 


|Mj(mAr)| = 


i/[Mf(mAr)] 2 + [Mf (mAr)] 2 + [Mf (mAr)] 2 


(7.37) 



Matrix multiplication by exp(— ArA(m)) is most easily done in real space. We have imple- 
mented the matrix multiplication by exp(— ArT) by means of the Fast Fourier transform 
instead of the usual checkerboard decomposition. In the discretized imaginary time version, 

T T exp ( - J h(s)ds) « exp(-Ar/i(ArL t )) . . . exp(-Ar/i(Ar)) 

= / ^(OIC > exp ( ]T ^[5^-! . . . B 1 ] iaJa ,$ a ,) < £1 (7-38) 

iajcr' 

Note that we choose to approximate the imaginary time integral of a function in an 
interval of length Ar by Ar times the value of the function at the upper limit of the 
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interval. The calculation of the normalization factor Z, Eq. (7.19), then becomes a simple 
exercise in Grassmann integration 

Z = J dnirfidnig) < -r]\i > exp ( - Y^™ + Y ^U^l^-i • • • ^iWv^V) < Z\v > 



= Det[l + Bi Bi... BjJ = Det[l + B> B{ . . . B 1 ^} (7.39) 

where B^ is the transpose of the matrix B m and in the last line we have used the periodicity 
of the exchange and scalar fields which implies B Lt = B and multiplied the expression by 
1 = Det[BT B r\ 

The remainder of the calculation of the Green's function proceeds along similar lines. 
Consider the case m > m'. One has 

T T exp ( - j 13 h(s)ds) « exp(-Ar/i(ArL t )) . . . exp(-Ar/i(Ar(m + 1))) 

= J dn(0\£ > exp ( Y U[B Lt B Lt ^ . . . B m+1 ] ia , ja ,C j( r) < f I (7-40) 



layer' 



with similar expressions for T T exp ( — h(s)ds^j and T T exp ( — J r h(s)ds^j. After straight- 
forward algebraic manipulations, one arrives at: 

1 t 

G(iam,ja'm') = -— J dp(i)dp(r])dp(p)dp(4>) exp ( - Y&M*) ex P ( ~ EOfcr) 





exp ( Y ViaPia) ex P ( 


Y PiMo\ 




1j*>) 




ia 


ia 


icr J, a' 






exp ( - Y Pi° M ljo 


<P*jo) exp 


( - Y ^Mi, 3 A%) 


(7.41) 




iaja' 




iajcr' 




M v 


= B Lt B Lt _i . . . B m+ i 






(7.42) 


M p 


= B m B m _i . . . B m / +1 






(7.43) 


M 4, 


= B m >B m '-i . . . B\ 






(7.44) 



The integral over the quadratic form is easily computed. The calculation in the case m <m! 
proceeds along similar lines and we just quote the final answer. The single particle Green's 
function is given by the following matrix expressions 



B T B T B T B T 1 4- B T B 1 B 1 B' 1 



T UT dT dT 



-1 



if m > w! 



G(m,rri) = < 1 _ 1 J (7.45) 



rT rT rT 
n m+l n m+2 ■ ■ ■ n m> 



1 + B m , +1 . ■■B Lt B 1 ...B 
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if m < m' 



Note that these expressions differ from those of Ref. [40] since these authors are using 
the standard formulation of the path integral with normal ordered operators, instead of the 
anti-normal ordering used here. The numerical calculation of the single particle Green's 
function, Eq. (7.45), is complicated by the fact that the matrix products entering the above 
expressions are seriously ill-conditioned [41]. I have found that the Gram-Schmidt matrix 
factorization algorithm of Ref. [41] can be unstable at sufficiently low temperatures and 
for strongly fluctuating molecular fields. I therefore use an alternative way to compute the 
Green's function which I have found to be particularly stable. It is based on the matrix 
identities: 

(1 + A 1 A 2 )- 1 = T 2 [l-T 1 -T 2 + 2T 1 T 2 ]" 1 T 1 (7.46) 
A 2 (l + A 1 A 2 y 1 = [l - T 2 ] [l - Tx - T 2 + 2TxT 2 ] "Vi (7.47) 
T l = (l + A i )- 1 (7.48) 

In practice we calculate G(m, 0) and G(0, to/). The above matrix identities can be used 
recursively to calculate these Green's functions in a stable manner. The matrices Ai, A 2 are 
taken as products of B T matrices. The number of B T matrices in a "block" is chosen to be 
the maximum number such that (1 + Ai), (1 + A 2 ) can be inverted without significant loss 
of precision. One can then obtain the Green's functions G(m, 0) and G(0, to') at a subset 
of values to, to'. The Green's function for the other values of to, to' can be obtained by 
propagating forward or backward in time. For example, G(m + 1,0) = [-B^ +1 ] _1 G(to, 0). 
We leave it to the reader to fill in the details. 
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FIGURES 
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FIG. 1. Diagrammatic expansion for the single particle Green's function. The dashed line 
connected at one end only to the solid line represents the interaction of the fermions with the 
dynamical molecular field and the brackets the average over the Gaussian distribution of these 
fields. The averaging over the distribution of molecular fields is carried out by pairing the dashed 
lines in all possible ways, each pairing giving a factor proportional to the two-point correlation 
function of the molecular field, the dynamical susceptibility x(q, iv n ) according to Eqs. (2.11,2.12). 
The lower part of the figure shows the pairings one obtains up to two spin or charge-fluctuation 
exchanges. 
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Eliashberg 
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Eliashberg ; k = (71,71/8) 
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Eliashberg ; k = (37r/8,37i/8) 
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FIG. 2. The one-loop approximation for the quasiparticle properties for both commensurate 
charge- fluctuations and antiferromagnetic spin- fluctuations. The tunneling density of states N(uj) 
is shown in (a) while (b) and (c) show the quasiparticle spectral function A(k, uj) for momenta just 
below the Fermi level, (b) shows A(k,uj) for a wavevector close to the Van Hove singularity and 
(c) shows A(k, uj) for a wavevector along the diagonal of the Brillouin zone. 



37 



Charge-fluctuations 
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Charge-fluctuations ; k = (71,71/8) 
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Charge-fluctuations ; k = (37r/8,37i/8) 




FIG. 3. Quasiparticle properties for a coupling to a scalar dynamical molecular field <£. The 
tunneling density of states N(oo) is shown in (a) while (b) and (c) show the quasiparticle spectral 
function A(k, to) for momenta just below the Fermi level, (b) shows A(k, to) for a wavevector close 
to the Van Hove singularity and (c) shows A(k, uo) for a wavevector along the diagonal of the 
Brillouin zone. 
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AF Spin-fluctuations 
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AF Spin-fluctuations ; k = (71,71/8) 
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AF Spin-fluctuations ; k = (37r/8,37i/8) 




FIG. 4. Quasiparticle properties for a coupling to an exchange vector dynamical molecular field 
M with antiferromagnetic correlations. The tunneling density of states N(u) is shown in (a) while 
(b) and (c) show the quasiparticle spectral function A(k, uj) for momenta just below the Fermi 
level, (b) shows A(k,uj) for a wavevector close to the Van Hove singularity and (c) shows ^4(k, uj) 
for a wavevector along the diagonal of the Brillouin zone. 
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k = (7C,7C/8) 
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k = (37i/8,3ti/8) 
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FIG. 5. Comparison of the one-loop self-consistent versus the non-perturbative calculation of 
charge- and spin-fluctuation exchanges, (a) shows the spectral function A(k, uo) for a wavevector 
close to the Van-Hove singularity and (b) shows A(k, uj) for a wavevector along the diagonal of the 
Brillouin zone. 
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Eliashberg 
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Eliashberg ; k = (71,71/8) 
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Eliashberg ; k = (37r/8,37i/8) 




FIG. 6. The one-loop approximation for the quasiparticle properties for ferromagnetic 
spin- fluctuations. The tunneling density of states N(oj) is shown in (a) while (b) and (c) show 
the quasiparticle spectral function A(k, to) for momenta just below the Fermi level, (b) shows 
A(k,uj) for a wavevector close to the Van Hove singularity and (c) shows A(k, uj) for a wavevector 
along the diagonal of the Brillouin zone. 
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FE Spin-fluctuations 
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FE Spin-fluctuations ; k = (71,71/8) 
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FE Spin-fluctuations ; k = (37r/8,37i/8) 




FIG. 7. Quasiparticle properties for a coupling to an exchange vector dynamical molecular field 
M with ferromagnetic correlations. The tunneling density of states N(u) is shown in (a) while (b) 
and (c) show the quasiparticle spectral function A(k, u) for momenta just below the Fermi level, 
(b) shows A(k,uj) for a wavevector close to the Van Hove singularity and (c) shows A(k,uj) for a 
wavevector along the diagonal of the Brillouin zone. 
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FIG. 8. First order vertex correction to the one- loop self-energy. The way spin is carried 
through the diagram depends on the kind of molecular field. A Pauli matrix is associated with each 
vertex in the case of a coupling of quasiparticles to magnetic fluctuations. In the case of exchange 
of charge fluctuations, each vertex simply carries a unit matrix, namely the spin orientation is 
unchanged at each vertex of the diagram. 
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k = (37i/8,3ti/8) 




co/t 

FIG. 9. Comparison of the single antiferromagnetic spin-fluctuation exchange approximation 
for the quasiparticle spectral function with a renormalized correlation wavevector n 2 e ^ = 0.35 
and the non-perturbative result with k 2 = 2. In both cases the dimensionless coupling constant 
9 2 Xo/t = 2. (a) and (b) show the quasiparticle spectral function A(k,Lo) for momenta just below 
the Fermi level, (a) shows A(k, oj) for a wavevector close to the Van Hove singularity and (b) shows 
A(k, to) for a wavevector along the diagonal of the Brillouin zone. 
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FE Spin-fluctuations 
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FIG. 10. Approximate quasiparticle mean-free paths /k for k = (tt,tt/8) and k = (3ir/8, 3ir/8) 
obtained from Lorentzian fits to the numerical results for the spectral function, (a) shows for 
coupling to charge- fluctuations, (b) to antiferromagnetic spin-fluctuations and (c) to ferromagnetic 
spin- fluctuations . 
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